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Abstract 

We show how Feynman diagrams may be evaluated to take advantage of recent 
developments in the application of Cutkosky rules to the calculation of one-loop am- 
plitudes. A sample calculation of gg — > gH, previously calculated by Ellis et al., is 
presented illustrating this equivalence. This example demonstrates the use of cutting 
rules for massive amplitudes. 
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1 Introduction 



One-loop calculations are important in the quest for discovering new physics beyond the 
Standard Model, especially for determining QCD backgrounds. A number of significant im- 
provements in the calculation of one-loop amplitudes, especially ones with massless particles, 
have been made in recent years. These methods, which have recently been reviewed in ref . HJ , 
include the use of helicity 0, color decompositions ||, string-inspired ideas recursion 
relations || ||, unitarity |7|, || |[ |10| and factorization [|TTJ. Besides allowing for the com- 
putation of certain infinite sequences of amplitudes, these techniques have been used in the 
computation of non-trivial amplitudes relevant for phenomenological applications [|1^, p~3|j . 

A technique that has proven to be especially useful is that of using unitarity to determine 
the functional form of the amplitudes. By working to all orders in the dimensional regular- 
ization parameter, e = (4 — D)/2, one can determine the coefficients of all integral functions 
that may appear [TJ|] . For certain amplitudes that satisfy a power counting criterion, such as 
supersymmetric ones, one can show that cuts with four-dimensional momenta are sufficient 
to determine the amplitudes through O(e ) || [J. The case of an amplitude with massive 
internal legs has been considered in ref. flOf . 

The efficiency of the unitarity method stems from the use of gauge invariant tree am- 
plitudes in performing calculations. With the use of helicity and color decompositions this 
can lead to enormous reductions in the complexity of intermediate expressions. This may 
be contrasted against calculations of individual Feynman diagrams that are gauge variant, 
where large cancelations must take place in order obtain simpler gauge invariant results. 

In this paper we will show how to evaluate Feynman diagrams so as to obtain efficiency 
of computation inherent in the unitarity technique. This is done by performing Feynman 
diagram calculations in a particular way so as to mimic the unitarity approach advocated in 
ref. [0. One may therefore start with Feynman diagrams, yet obtain the advantages of the 
Cutkosky approach. This is useful since it explains the unitarity method from a conventional 
Feynman diagram viewpoint. This is of interest in, for example, the extension of unitarity 



methods beyond one loop, such as in ref. [15 



As an explicit example we have calculated the one-loop helicity amplitudes for gg — > gH . 
This example has the desirable features that it contains both colored and colorless lines and 
massive internal and external legs. This example is useful since it illustrates features beyond 
purely massless calculations, where most attention has been concentrated M. Also it has 



already been evaluated via conventional methods by R.K. Ellis et al. [|16[], allowing for an 



easy comparison to the correct result to be made. We find agreement with their result. 



2 Review of basic tools 

We now briefly review some of the techniques that we shall use in this paper. Two important 
techniques that we utilize are spinor helicity and color decomposition. The reader may 



consult review articles [171 for further details 



1 



2.1 Spinor helicity 



In explicit calculations with external gluons, it is usually convenient to use a spinor helicity 
basis which rewrites all polarization vectors in terms of massless Weyl spinors \k ± ). In 
the formulation of Xu, Zhang and Chang the plus and minus helicity polarization vectors 
are expressed as 

where k is the momentum of the gluon and q is an arbitrary null 'reference momentum' 
which drops out of the final gauge invariant amplitudes. We use the convenient notation 

(kr\k+) = (ij) , (kt\kj) = [ij\. (2.2) 

These spinor products are anti-symmetric and satisfy (ij) [j i] = 2ki ■ kj. One useful identity 



is 



(a\ fyk; q) \b) = ± 7^rM [I^X^I + I^WIJ \ b ) > (2-3) 

where either (a\ or |6) are spinors with four- dimensional momenta. 

It is convenient to choose a dimensional regularization scheme that is compatible with the 
spinor helicity formalism. We use the four-dimensional helicity scheme || which is equivalent 
to a helicity form of Siegel's dimensional reduction scheme ||18|| . The conversion between the 
various dimensional regularization schemes has been given in ref. [|], [HJ . 



2.2 Color decomposition 



Color decompositions M have been extensively discussed in review articles |T7[. Here we 
only present the color decomposition for the gg — > gH amplitude. 

For this one-loop amplitude the decomposition is 

^i-ioop = ^|^^ TMl)TM2)T a CT(3) ] A{{a{l),cr{2),a{3),H), (2.4) 

(7 

where we have taken the loop to be in the fundamental representation. The sum over a 
includes all cyclic permutations of the indices a{n) (i.e. a G {(1,2,3), (1,3,2)}) and the T a 
are fundamental representation color matrices (normalized so that tr (T a T h ) = 5 ab ). We have 
abbreviated the dependence of A{ on the outgoing momenta kj and helicities Xj by writing 
the label j alone. We have also explicitly extracted the coupling and a factor of /i^ e from A{, 
where fin is the renormalization scale. For adjoint representation loops there is an analogous 
decomposition with up to two color traces in each term. 

It is convenient to further decompose the partial amplitude in terms of primitive ampli- 
tudes, 

A{(a(l),a(2),a(3),H) 

= A 4 {a{l), (7(2), (7(3), H) + At(<7(3), (7(1), (7(2), H) + A 4 (a(2), (7(3), (7(1), H) . 

(2.5) 
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In the primitive amplitudes appearing on the right-hand-side we may treat the Higgs as car- 
rying color charge (in the adjoint representation), which couples to gluons; it is not difficult 
to show in the permutation sum all such additional terms cancel. This is convenient since 
these primitive amplitudes are gauge invariant and have a fixed cyclic ordering of external 
legs. A similar decomposition in terms of 'primitive amplitudes' can also be performed for 
the case where some of the external particles are in the fundamental representation p0| . 



3 Feynman diagrams 

Feynman diagrams, which have been essential for evaluating terms in the perturbative se- 
ries for scattering processes, suffer from a number of well known computational difficulties. 
Firstly, individual Feynman diagrams are inherently gauge dependent; gauge invariance is 
restored only after summing over all diagrams. Secondly, when evaluating loop momentum 
integrals using traditional techniques, such as Passarino-Veltman reduction , one encoun- 
ters a large number of spurious denominators. These spurious denominators typically involve 
Gram determinants, det(fcj • kj) raised to a power. Enormous cancelations must then take 
place before one obtains relatively simple results. 

When using conventional techniques to evaluate Feynman diagrams one reduces all tensor 
integrals into sums of scalar integrals multiplied by reduction coefficients. (By a tensor 
integral we mean loop momentum integrals with powers of momentum in the numerator and 
by a scalar integral we mean integrals with no powers of loop momentum in the numerator, 
see appendix With the knowledge of how this reduction works, given any one-loop 
amplitude we can list all the scalar integrals that could potentially enter. Thus we have 

Amplitude = V] Feynman diagrams — > CjJj , (3.1) 

i 

where the subscript i is some index that labels the scalar integrals that contribute to this 
amplitude. The heart of the calculation is in evaluating all the coefficients of the scalar 
integrals, i.e. the c^'s. 

A conventional Feynman diagram technique like Passarino-Veltman (PV) reduction eval- 
uates all the necessary coefficients simultaneously. If one could obtain these coefficients in a 
more efficient manner one would greatly improve the ability to calculate. 

Some of the features that are useful in improved calculational methods are: (1) that the 
building blocks be gauge invariant combinations of diagrams instead of individual diagrams; 
(2) that wherever possible tensor integrals should not generate large powers of Gram deter- 
minants in denominators; (3) that large cancelations between diagrams be avoided as much 
as possible; (4) and that calculations can be recycled when performing new ones. 
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Figure 1: A cut of a general n-point one- loop amplitude. 



3.1 Review of the unitarity method 

The cutting method reviewed in ref. |l| has these above properties, 
writes down the cuts in a given channel, depicted in fig. [I], 



In this method one 
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(3.2) 

This equation is valid only for the cut channel under consideration. The loop momentum is 
p; l\ and £2 are the momenta of the cut propagators, with masses mi and 1712 respectively 
and rii + 112 = n + 4. One then reconstructs the complete amplitude by finding a function 
which has the correct cuts in all channels. 

Since the loop integral in eq. ( pT2| ) is composed of products of on-shell tree amplitudes, 
the basic building blocks are gauge invariant. These tree amplitudes are inherently easier 
to calculate than loop amplitudes and in many cases have been previously calculated. A 
key feature of the cutting method is that a specific subset of scalar integral coefficients are 
obtained in any given channel. Then one repeats this process for various channels thus 
computing the coefficients for all the scalar integrals. There are integrals that do not appear 
in any channel, namely the tadpole and scalar bubble integrals, 



D=i-2e 



and 



lim /, 

fc 2 ->0 



D=4-2e 



(3.3) 



However, with prior knowledge of the infrared and ultraviolet behavior of an amplitude in 
many cases the coefficients of these integrals can be fixed, as discussed in ref. [ [10|j . 

It is generally convenient to express the tree amplitudes in eq. ( |3.2| ) using the spinor 

Then one can evaluate the right-hand side of 
a\ ■ ■ -l-i ■ ■ ■ (-2 ■ • ■ |&)) an d then commuting 
i ■ kj, which may be rewritten as 
differences of inverse propagators and terms which do not depend on the loop momentum. 
This type of reduction procedure has been previously used in a number of references || ||. 
PH[ P2] |. We shall make use of it in the next section. 



helicity formalism discussed in section ^1 

I2|) by constructing spinor strings of the form 



eq. flo.^i ; uy LAjiiDLiu^umg, opmui ouiiiig,D ui tuc luiui \u,| ■ • ■ p x 

the li towards each other. This generates terms of the form 
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Figure 2: The contribution from a vector in the loop is separated into the difference of a 
vector and scalar plus a scalar so that the difference satisfies the power-counting criterion. 



When using spinor helicity one must take into account that the momenta of the inter- 
mediate on-shell lines are in (4 — 2e) dimensions while the helicity formalism is defined in 
four dimensions. If one drops the (— 2e)-dimensional components of the momenta, one could 
make an error in rational functions arising from e/e terms. For five- and higher-point cal- 
culations it can be more convenient to use four- dimensional momenta in the cuts, but then 
to reconstruct the missing rational functions using collinear and multi-particle factorization 



properties [II], |TJ. For four-point calculations it is convenient to keep terms to all orders in e, 
since factorization is insufficient to fix all rational functions. 

The cases of fermion or scalar lines with (4 — 2e) momenta in the helicity formalism 
have been extensively discussed in ref. and do not present any difficulties. Although 
it is unnecessary for the specific calculation gg — > gH, a convenient strategy for handling 
channels with gluon lines is to subtract and add equivalent diagrams with internal gluon lines 
replaced by scalar or fermion lines, as depicted in fig. |2|. In some cases one uses fermions so 
as to group particles together with supersymmetric partners. Such combinations of diagrams 
satisfy the power counting criterion that n-point tensor integrals should have (for n > 2) 
at most n — 2 powers of loop momentum in the numerator of the integrand; two-point 
integrals should have at most one power of loop momentum. The required rearrangement 
of the diagrams to make the power count apparent is discussed in ref. 0. In particular, for 
the example portrayed in fig. ^| one should use background field Feynman gauge . The 
origin of this power counting criterion is from the requirement that all contributions due to 
the 0(e) parts of loop momenta should not interfere with ultraviolet poles to yield O(e ) 
contributions. (One can also show that infrared divergences never interfere with the 0(e) 
parts of loop momenta |§.) This criterion also holds for the general massive amplitudes 
discussed in this paper. This is useful since one can then take the sewed gluon lines to have 
four-dimensional momenta without any errors through O(e ). 



3.2 From Feynman diagrams to the unitarity method 

We now describe how to perform a Feynman diagram calculation so as to mimic the cutting 
method. In the next section we give an explicit example involving both internal and external 
masses. The basic idea is rather simple. 

Given that we know all the scalar integrals that could possibly enter a one-loop amplitude, 
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suppose we target those scalar integrals that have a pair of scalar propagators 



i/(£i-ml) and i/(t 2 - m 2 2 ) , (3.4) 

where i\ is the momentum of the scalar propagator between external legs i\ — 1 and i\, and 
£2 is between i 2 and i 2 + 1 of a general n-point amplitude, where the indicated propagators 
are each of a specific particle type. In order to generate a scalar integral with a specific prop- 
agator one must necessarily have evaluated a Feynman diagram that has the corresponding 
propagator. 

We can define a channel (this channel will correspond to the cut-channel when the cut 
is performed on the indicated legs) as a subset of the terms in the full amplitude; the terms 
being the scalar integrals which include the pair of scalar propagators mentioned above 
together with their corresponding coefficients, the q's. 

In each Feynman integral (the tensor integral corresponding to each Feynman diagram) 
we explicitly extract a denominator that includes the pair of scalar propagators under consid- 
eration. Then if we impose the on-shell conditions i\ — m\ = l\ — m\ = to the numerators 
of all the Feynman integrals contributing to the amplitude then the result would yield ex- 
actly what we defined as a channel in the previous paragraph. The potential errors when 
imposing these on-shell conditions would necessarily be of the form 

, (3.5) 



where n > 1 and/or m > 1, and D = 4 — 2e. Since (l\ — mf) and (£ 2 — vn^) are themselves 
inverse scalar propagators, anything of the form of eq. ( |3.5| ) could only generate terms which 
are not in the channel under consideration. This is because one or both of the inverse 
propagators would cancel a scalar propagator in the denominator and generate a scalar 
integral not included in the specified channel, i.e. that does not include both of the indicated 
scalar propagators. Setting i\ — m\ = i\ — m\ = only affects the calculation of those terms 
not in the specified channel. Thus the evaluation of the Cj's that are in the specified channel 
are not affected by imposing the on-shell condition. 

By systematically stepping through a sufficient number of channels one can fix all the 



coefficients, the Cj's in eq. ( |3.1|) . In this way we can reconstruct the full amplitude. For some 
coefficients the corresponding integrals appear in multiple channels; for consistency these 
must coincide. 

So how does using the on-shell conditions i\ — m\ = £ 2 — m 2 = in the numerator help 
us in evaluating the Feynman diagrams? All Feynman diagrams in a specific channel can 
now be 'factorized' using their on-shell legs into the sewing of two Feynman tree diagrams 
(see fig. §). 

In figure fig. |3] the 'factorized' Feynman graph only contains intermediate on-shell fermions 
for the gg — > gH amplitude. But in general this 'factorization' process can be performed 
on any one-loop Feynman diagram. With on-shell scalars or fermions it is clear that this 
can be done. With intermediate massless gauge particles one might worry that one has 
forgotten about the contributions of diagrams with ghosts on the right hand side of fig. 0. 



6 



H 













mnnnr 









H 

_y U D i 

helicity 

J (271)° 



X ■ 



(£?-m?) (el-mi) 



■ x 



t]-ml = 



Figure 3: An example of a Feynman graph for the gg — > (7// amplitude that contributes to 
the s-channel, 'factorized' into the sewing of two trees. The dashed lines on the left hand 
side indicate the on-shell propagators and the helicity sum is over the polarizations of the 
'factorized' legs. 



We know, however, that the net effect of the one-loop diagrams with ghosts is to cancel the 
unphysical polarizations of graphs with corresponding gauge particles. This cancelation can 
be made apparent using the background field Feynman gauge ||23|| . However, once one has 
'factorized' a Feynman graph on a pair of propagators that includes a gauge particle then 
the contribution of the ghosts is not necessary as the sewing of the two on-shell Feynman 
trees manifestly has no unphysical gauge degrees of freedom. 

If one sums over all the Feynman diagrams that contribute to a specific channel and we 
impose the on-shell condition then the equivalent equation to fig. |3| would have on-shell tree 
amplitudes sewed together instead of on-shell Feynman trees: 

A,(M,...,n)i channel = £ f^:uv., !2 



xA^K,Ul,.,«i-l, 



-5*1, 



t\-m\=l\-m\=§ 

(3.6) 

The sum on the right-hand side runs over all helicity states that can propagate across 
the 'cut' legs. Using the above equation one can compute all the coefficients of the scalar 
integrals that appear in this channel. 

The above equation, eq. (|3.6|) , is the unitarity method mentioned earlier in this section. 
Originally the method was inspired by Cutkosky rules, however, from the arguments pre- 
sented here it is clear that one can evaluate Feynman diagrams in a sophisticated manner so 
as to take advantage of the benefits inherent in the unitarity method. Instead of performing 
the calculation of the scalar integral coefficients in parallel, we perform the calculation by 
targeting simultaneously only those integrals in a given channel. 

The notion of channels presented in this section is more general than that of just consid- 
ering Cutkosky rules because we may use the on-shell conditions on any propagator which 
does not collapse in the integral function under consideration. (Indeed, this observation has 
already been used in calculations of Z —>■ 4 partons ||13|| .) In general we could construct 
a channel where we use the on-shell condition for any number of propagators (it is even 
possible to construct a channel where the on-shell condition is used for a single propagator); 
however, it is normally convenient to use only the two-propagator channel construction. 
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4 Calculation of gg —> gH 



In this section we apply the procedure outlined in the previous section to the specific example 
of gg — > gH. This amplitude has one massive external leg and all the Feynman diagrams 
have a massive quark running in the loop. See fig. ^ for the color-ordered Feynman diagrams 
contributing to this amplitude. 




Figure 4: The color-ordered Feynman diagrams contributing to the amplitude gg — > gH. 

In fig. |] there are two Feynman diagrams that involve a Higgs-Higgs-gluon coupling which 
does not exist in the Standard Model. As mentioned in section \2.2[ although these diagrams 
cancel from the final result, they make the primitive amplitudes (which have a fixed ordering 
of external legs) gauge invariant. Of course, one does not need to include the extra diagrams, 
but then one would need to sum over the permutations with all possible external orderings 



of the Higgs, as in eq. Q2.5Q , to recover gauge invariance. 

As mentioned in section [|by observing the Feynman diagrams in fig. f|we can write down 
all the scalar integrals that can enter the amplitude, see fig. The integrals are two-, three- 
and four-point scalar integrals that have massive internal propagators and some combination 
of massive external kinematics. One-point tadpole integrals do not appear since they are 
forbidden by power counting. (When using Passarino-Veltman reduction one could encounter 
such integrals, but they must cancel in the final expressions for this amplitude since one does 
not encounter them when using Feynman parameterization.) The scalar integrals are found 
in appendix [A.l. 



In addition to the use of color decomposition, section pl2| , it is also convenient to decom- 
pose the primitive amplitudes into all possible external helicity configurations. So the objects 
that need to be computed are the color ordered gauge invariant primitive amplitudes. There 
are three such objects; they are A 4 (l + , 2+ 3 + , H), v4 4 (l + , 2 + , 3", H) and A 4 (l + , 2~, 3+, H). 
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Figure 5: The — > gif amplitude parameterized in terms of all the scalar integrals that 
might enter. The external legs are labeled clockwise starting at the top right. 

The other primitive amplitudes (e.g. Ai(l~, 2~, 3 + , H)) can be obtained from the above 
three primitive amplitudes using parity and relabelings. 



Since this amplitude has already been calculated by Ellis et al. [16j the purpose of the 
calculation is to explicitly demonstrate that we may obtain the correct results by applying the 
unitarity method to massive amplitudes, starting from Feynman diagrams. For illustrative 
purposes, only the calculation of A 4 (l + , 2 + , 3 + , H) is explicitly presented here - the results 
of the calculation of the other two primitive amplitudes are shown for completeness. 



4.1 The Al(1+ 2+ 3+, H) primitive amplitude 

In order to compute the primitive amplitude Ai(l + , 2 + , 3 + , H) we will construct it from 
multiple channels. We will see that it is sufficient to only compute the s- and t-channels. 
The Mandelstam variables are defined as s = (k± + A^) 2 and t = (k 2 + k^) 2 . These channels 
can be seen in fig. || 

We begin by first constructing the s-channel by sewing together the two tree amplitudes 
on the left- and right-hand side of fig. ^|a. This corresponds to targeting the coefficients, ci, 
C2, C3 and C6 in fig. |]. The tree amplitude on the right of fig. |]a is given by |10 



where ui± = |(1 ± 75) is the helicity projection operator and the tree amplitude on the left 
is given by 

AT-i-L* 3+, H, LO = -™- ( ffil - y - 2 Ll) + 2 -^4(-L 2 \ - L l} ) . (4.2) 



V2v \ (e 1 + h) 2 - ii 



s — m 



Here v = is the vacuum expectation value of the scalar Higgs field ]2"1 ], m is the mass of 
the fermion in the loop and m# is the scalar Higgs mass. Capitalized momenta are (4 — 2e)- 
dimensional, lower case momenta are four- dimensional and /i is the (— 2e)-dimensional part 



9 



H 1 




Figure 6: The s- and t-channels of the gg — > gH amplitude. The on-shell lines are fermions 
since at tree level quarks are the only particles to couple both to gluons and to the scalar 
Higgs boson. 



of the loop momentum. (See appendix [A] and appendix 0.) In these tree amplitudes we 
follow the conventions of ref. [10] for massive spinors, summarized in appendix [B]. In this 
convention the distinction between fermion and anti-fermion spinors is suppressed. 

As depicted in fig. |]a, for the s-channel the momenta of the on-shell legs are 



P, 



p — ki — k 2 



(4.3) 



The second tree amplitude, eq. (|4.2| ), involving the Higgs has been constructed from the tree 
level Feynman diagrams including the additional fictitious Higgs-Higgs-gluon coupling (this 
is the standard scalar-scalar- vector coupling of scalar QED ]f24| ) mentioned earlier. 

So sewing the two tree amplitudes together gives us the s-channel of the amplitude 



V2mk 4 -et [12] r d^ 2e P M 



s— channel 



V(S — m 



Hj 



(12)7 (2k) 
m [12] 



4-2e 



V2v(l2)J (2 



V 3 

d 4-2ep 



i4-2e 



(4.4) 



where the denominators are 



fj 2 — m 2 



2 2 

jj, — m 



and the numerators are 



i — k\) 2 — fj 2 — m 2 
i — ki) 2 — /i 2 — m 2 



i — kx — k 2 ) 2 - /i 2 - m 2 
i — kx- k 2 ) 2 - /i 2 - m 2 



i + k^) 2 — /i 2 — m z 



M = (-AK^ + m)! -L 2 ){-L 2 \ -L^ 

M 2 = (-L 1 \u+(i/L + m)\-L 2 )(-L 2 \t+(t 4 -2m)\ - L x ) . 



(4.5) 



(4.6) 
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The two terms in eq. ( |4.4| ) correspond to the two terms in eq. ( ^4.2| ). 

First, concentrating on the numerator Mi we can write it in the form of a trace 

M = tr+[Ql + m)(/ 2 + !/l - m)(A + - m)] , (4.7) 

where we use the notation tr±[. . .] = tr[u;±...]. When expanding out the trace we first 
concentrate on the rt's since the final expression for the numerator can only depend on fi 2 
since the (— 2e)-dimensional integral over /i vanishes for a integrand with a odd number of 
/i's. Thus, 

M = tr + [(m/ 2 - fi 2 - m 2 ){t x - m)] 

= 2m(m 2 + fi 2 ) + 2m£ 1 - t 2 (4.8) 
= 4m(m 2 + fi 2 ) — ms . 

The numerator M 2 can also be written in the form of a trace 

M 2 = tr + [fyt + m)(/ 2 + fjt - m)^(h - 2m) (A + ^ - m)] . (4.9) 

So again first concentrating on the it's we can write 

Af 2 = tr + [(^ 2 + m/ 2 - /i 2 - m 2 )4(# 4 - 2m) (A + ^ - m)] 

= tr + [(m/ 2 - /i 2 - m 2 )4(# 4 - 2m) - m)] - fi 2 tr + [/ 2 ^(# 4 - 2m)] (4.10) 
= 2m(m 2 + /i 2 )(24 + 2i 2 + fc 4 ) • 4 + mtr+^+^J . 

Recalling that l 2 = i\ — k\ — k 2 and choosing the reference momenta q% = k\ (so that 
fci-4 = 0) then 

M 2 = 8m(m 2 + fi 2 )i x ■ 4 + 6m(m 2 + fi 2 )h • 4 + mti+fafi&Ji] • (4.11) 
We can evaluate the remaining trace in the expression for M 2 by using eq. 



tr+t/,4^] = tv^UM = ^!|<:r|444|3-} , (4.12) 

where for brevity we represent external momenta by their indices and neglect slashing the 
momenta inside the inner-product. Using / x / 2 = /i(/i — fi — $2) = m 2 + /i 2 — /i($i + fz 2 ) 
and momentum conservation the trace can be expressed as 

tr+^4^] = ^| {(m 2 + /, 2 )(1-|4|3-) + (r|4^4|3-)} . (4.13) 

But since 



, (1~|4|3-) 

k --^--JiU' (4 - 14) 



M 2 can be expressed as 



M 2 = 8m(m 2 + fi 2 ){h • 4 + h • 4) + v^m ^ 1 ^ |3 } . (4.15) 

(13) 
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The integral comprised of the numerator J\f 2 now only involves at most a single power of 
loop momentum i in the integrand. This can be readily evaluated with the use of Feynman 
parameterization. The Feynman parameter shift is 



?i = q + a 2 k x + a 3 (kx + k 2 ) - a 4 k 4 



(4.16) 



where q is the shifted loop momentum. After Feynman parameterizing the integral the terms 
in M 2 involving l\ can be evaluated by performing the momentum shift and next dropping 
terms containing a single power of q which vanish. So the two terms in M 2 involving l\ 
become 

(l~\q + a 2 ki + a 3 (fci + k 2 ) - a 4 fc 4 |3~) 



^(13) 



— (a 3 + a 4 )k 4 .e^ , 



and 



(1-|4^4|3~) = a 2 <l-|414|3-> + a 3 <l-|4(l + 2)4|3~) - a 4 (l- 1444|3~> 
= [(t — m 2 H )a 2 — sa 3 — m^a4](l~|4|3 - } , 



(4.17) 



(4.18) 



after performing the momentum shift. 

As defined in appendix [A.3| , the integral function I n [aj\ is a Feynman parameterized form 
of the loop integral with a factor of in the numerator. Thus after Feynman parameteriza- 
tion eq. (4.4) becomes 



A 4 (l + ,2 + ,3 + ,H) 



s~ channel 



— I 



V2m 2 k 4 ■ ef [12] f 
v{4tt) 2 -^ (12) { 



m 



2~ J 3 
H 



4 4) [4(m 2 + fi 2 ) -s]-{t- m 2 H )I 4 [a 2 ] 



+ sl 4 [a 3 ] + m 2 H I 4 [a 4 ] - AI 4 [(m 2 + /i 2 )(l - a 3 — a 4 ) 



where I, 



D=4-2er,,2 



[[i 2 ai 



-el. 



D=e-2e\ 



and J 3 is defined in appendix [A]. The dimension label 



s— channel 

(4.19) 



on (4 — 2e)-dimensional integrals has been suppressed. Then using the recursive formula for 
In[o>i] m eq. (|A.14 ) and eq. (4.14) one can write the s-channel of the primitive amplitude as 



A 4 {l + ,2+,3 + ,H) 



m 



st 



t — u 



4 4) [4K + /i 2 



s— channel 



(4vr) 2 - e v (12) (23) (31) [2t(s-m 



m 



m 



H r(2) 



1 



2st 



I^>[A(m 2 + /i 2 ) - m 2 H ] - -I 4 [A(m 2 + // 



m 



Hi 



) s— channel 

(4.20) 

where u = m 2 H — s — t, also noting that J 3 2 ' ) is a scalar triangle integral dependent on the 
kinematic invariant s and m 2 H , while J 3 4 ^ is a function of s only. Then using eq. (|A.11|) the 
s-channel becomes 



A 4 (l + ,2+,3 + ,#) 



s-m H (2 ) 



m 



st 



s-channel (Aw) 2 ^ V (12) (23) (31) 



(4m 2 — m 2 H ) 



t — u 



2t(s - m 2 H ) 



2st 



-4e 



t — u 



2t(s — m 2 H 



r (4)D=6-2e 
2 \ 1 3 



m 



H r (2)D=6-2e 



2st 



T D=6-2e 
1 4 



s— channel 

(4.21) 
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Now that we have the s-channel of this amplitude we need to calculate the i-channel 
(corresponding to targeting Ci, C4, C5 and C7 in fig. |J) in order to reconstruct the full primitive 
amplitude. However due to the symmetry of the all-plus primitive amplitude it is not 
necessary to explicitly perform this calculation as the s- and t-channels are related up to an 
overall sign by a simple exchange of legs 1 and 3 which is equivalent to exchanging s and t 
in eq. ( |4.21| ) (and also switching the 1 and 3 in the inner products in the pre-factor of the 
expression). Thus the t-channel of A 4 (l + , 2 + , 3 + , H) is given by 



A 4 (l + ,2+,3 + ,#) 



m 



st 



t— channel 



(4tt) 



(4m 2 — m 2 H ) 



s — u 



2s(t-m 2 H ) 



rn 



2st 



h AS) 

1 3 



2 



4e 



s — u 



2s(t 



m 



H 



v (12)(23)(31) 

_j(l)D=6-2e _ t - m 2 H (3)D =6 -2£ _ \jD= 

) 3 2st 3 2 4 



-2c 



(3) 

again recalling that 1% is a scalar triangle integral dependent on the kinematic invariant t 
and m 2 H) while I3 is a function of t only. It is also useful to observe that the coefficient of 
the I4 integral in the t-channel matches that of the s-channel as expected. If we refer back to 
fig. |5], where all the scalar integrals that could possibly enter were listed, we notice that the 
coefficients, eg of ^{mj^) and eg of -^(O) have not yet been fixed. However, there are no other 
UV infinities in either the s- or t-channel expressions (eqs. ( ^4.21|) and ( f4.22|) ), so therefore 
since this primitive amplitude is UV- finite (since it vanishes at tree level), it is not surprising 
that these coefficients are zero. This may be explicitly verified by direct construction of the 
channel corresponding to using the on-shell condition on the two propagators adjacent to 
the massive external leg, showing that both c§ and eg are individually zero. We have not 
evaluated this channel, since this example is only for illustrative purposes. 

Thus with the two channels of the primitive amplitude we can construct the full primitive 
amplitude A^(l + , 2 + , 3 + , H), which is 



t— channel 

(4.22) 



A, 



;i+,2 + ,3- 



H) 



m 



st 



(4yr) 2 - e v (12)(23)(31) 



(4m 2 



H) 



s — u 



r(l) 



2s(t - m 2 H Y 



s-m H (2) 
2st 3 



m 



H A3) 



2st 



+ 



t — u 



2t{s 



in 



H) 



r(4) 
3 

2 



-4e 



u 



r (l)D=6-2e 



s - m H A2)D=e-2e t - m H ( 3 )D=6-2e 
: — : :: — : — 



2st 



2st 



+ 



2s (t — m 2 h 
t — u 
2t(s-m 2 H ) 



r (4)D=6-2e 



2 4 



-2c 



(4.23) 

Note that this expression is valid to all orders in e with four-dimensional external momenta. 
Since this amplitude has no ultraviolet or quadratic divergences, the primitive amplitude is 
well defined in the limit e —>■ 0. Therefore eq. ( |4.23| ) can be simplified, since in this limit with 

I and e/ 4 D=6_2<: — > 0. Also with the intention of simplifying 



the use of eq. (Q, el? =6 ~ 2e 
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this expression we can use eq. ( |A.5| ) explicitly. Thus 

i m 2 st 



A 4 (l + ,2 + ,3 + ,H) 



(4tt) 2 v (12)(23)(31) 
+ ^I 3 ( m 2 H )-±I 4 (s,t,m 2 H 



(Am — m H ) 



his) 



s — m 
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2u(m A H — ts) 



t — m 2 H 



ts(s — m 2 H )(t — m 2 H ) 
where all the integrals in this equation are in four dimensions. 



h{t) 



(4.24) 



4.2 The remaining primitive amplitudes 

In a similar fashion to the all-plus helicity case the primitive amplitudes corresponding to 
the remaining two helicity configurations may be calculated. The results are 



A 4 (l + ,2 + ,3-,H) 



m 2 [12] 3 t 



Aum 



(An) 2 v [13] [23] s \(t-m 2 H ) 2 



h{t) - 



Aum 2 m 2 



H Um 2 



t(t - m 2 H ) 2 



H) 



+ (Am 2 - s) 
2u 

~ (t - m 2 H ) 2 



^h(s) 

s — m z H 



hit) - h(m 2 H/ 



t — m 2 H 



+ 



h(t) - 



ni 



t(t — m 2 H 



-h(m 2 H ) - -I 4 (s,t,m 2 H) 



2u(t - s) 



t(s-m 2 H )(t-m 2 H ) 



(4.25) 



and 



+ 



,3 + ,H)~- 
u — Am 2 



m 



H 



As 
+ — + 

u 



i m 2 [13] 3 st 
(4^)2 V [12] [23] ^ 
Aum 2 



2 4ts 

u — 12m 

u 



m 2 H ) 2 



I A (s,t,m 2 H 
Aum 2 



+ m H 



(u — Am 2 )(st — u 2 ) 

st(s — m 2 H )(t — m 2 H ) 

2(Zu + 2t) { ~ ~ a 

h(s) - I 2 (m 2 H/ 



_ , . \u — Am 2 At 

h(s) + — + — + . 2 s 2 

t — m H u (t — m u y 

Aum 2 Aum 2 



h(t) 



u 



nr 



2 H? 



s(s — m 2 H ) 2 t(t — m 2 H ) 2 
2(Zu + 2s) r ~ . , ~ , 2 , 



2u(u 2 - st) 



st(s — m 2 H )(t — m 2 H ) j 
(4.26) 

where the h's are bubble integrals with their 1/e poles removed (since the poles cancel in 
the difference between the two /2's) and all integrals are in four dimensions. 

Thus in equations ( [4.24|j4T2"5|14.26| ) we have the three necessary primitive amplitudes 
A 4 (l + ,2+,3 + ,iJ), A 4 (1 + 1 2 + ,3~,H) and A 4 (l+, 2~ 3+ H) which are needed for the full 
gg — ► g H amplitude. We have compared these amplitudes to those of Ellis et al. |Tj| and 
they agree. This comparison is performed by dotting the tensor expressions of Ellis et al. 
with spinor helicity polarization vectors. 
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5 Summary and discussion 



In this paper we have explained the relationship of the unitarity method reviewed in ref. 
to Feynman diagram calculations. One may calculate Feynman diagrams in a sophisticated 
manner so as to take advantage of the benefits of the unitarity method. We have illustrated 
the usefulness of this method for the case of gg — ► gH. This example is general in that it in- 
corporates both massive internal and external kinematics as well as a combination of colored 
and colorless external particles. It should be evident from the explicit calculation presented 
in this paper that this method lends itself to making complicated analytical calculations 
more tractable. This is because in this approach one is always dealing with gauge invariant 
constructions, which are generally significantly simpler than the gauge dependent Feynman 
graphs in a conventional calculation. In the unitarity method one targets the coefficients of 
a subset of scalar integrals contributing to the amplitude by computing a gauge invariant 
combination of Feynman diagrams. 

The notion of a channel as defined in this paper is more general than just the cut legs of a 
Cutkosky calculation, where one is nominally computing the discontinuity of an amplitude in 
a certain cut channel. The reason it is more general is that one defines a channel as a subset 
of the terms in the full amplitude; those terms being the scalar integrals that include some 
pre-defined combination of internal propagators. This notion of a channel has already been 
useful in the evaluation of Z —>■ 4 partons |0|. It was also helpful in two-loop calculations 



in N = 4 supersymmetric Yang-Mills [15 



One of the themes of this paper has been to explain the usefulness of the unitarity method 
to multi-point one-loop calculations. The ideas discussed here will potentially be beneficial 
to the extension of these types of calculations to multi-loops. Initial steps have already been 
made fl5] , and more attention is being devoted to these pursuits. 
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A Integrals 



In this appendix we collect expressions for the integrals [21], |26| we use in this paper. 



The D- dimensional loop integrals considered in this paper are defined by 

d Dp pai... pa m 



I"[p al • • • P ar 



- 1 (4tt 



.D/2 



(A.l) 



(27r) D (P 2 -m 2 )...((P-EL 1 

where P is a D-dimensional momentum and the external momenta, hi, are four- dimensional. 
Our convention is to suppress the dimension label on I n when we are dealing with aD = 4— 2e 
integral, unless otherwise indicated. We also suppress the square brackets when the argument 
is unity. Such integrals are referred to in the text as scalar integrals, while tensor integrals 
have m > 1. 
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Given a four-point amplitude, we have two conventions for labeling the one-, two- and 
three-point integrals. The first, and more familiar one is to explicitly give the kinematic 
invariant upon which the integral depends. For example, h(s) is the (4 — 2e)-dimensional 
bubble that has an invariant mass square of s = (ki + /c 2 ) 2 flowing through its external legs. 
The second convention for labeling the lower point integrals is to indicate, by a raised (i) 
that the internal propagator between legs i — 1 (mod 4) and i has been removed. If two 
indices are present as in then two propagators have been removed. This follows the 
labeling convention of refs. p6| , |T0] . 



A.l The D = 4 — 2e scalar integrals 

The bubble with an external kinematic invariant s is 



his) 



r(2,4) 



/ 2 (0) + 2 + xlog 



x — 1 

X + 1 



+ 0(e) 



(A.2) 



where x = yl — Am?/s. For external massless kinematics, this may be written in closed 
form as 



- 2e r(i + e) 



/ 2 (0) = m 

t 

The scalar triangle integral can either be a function of just one kinematic invariant s 



his) 



or of two invariants s and m 2 H 



r(4) 



X + l 
X — 1 



+ 0(e) 



h(s,m 2 H ) 



r(2) 



sl 3 (s) - m 2 H h(m 2 H ] 



s — m 



+ 0(e) 



(A.3) 



(A.4) 



(A.5) 



li 



While the box with a uniform internal mass and a single massive external leg, with invariant 
mass-square m 2 H is 



I 4 (s,t,m 2 H ) = - — 



um 2 m 2 



um 2 m 2 



where u = m 2 H — s — t, 



H(X,Y) = 



x + — X- 



St 



As 



St t 



st ' m 2 H 



+ 0(e), (A.6) 



+ Li 5 



y-x+ 



X -I 



y 



+ Li 2 ( 



x--y 

X-\- 

y-x. 



(A.7) 



with x± = |(1 ± VI -4A), y 



\/l — 4Y) and the dilogarithm is defined by 



Li 2 (a;) = - J a ln(l - xt)/t. 

The scalar integral expressions quoted above are only valid in certain analytic regions. 
These can be analytically continued, but are not presented here since related expressions^] 
are presented for all analytic regions in ref. [16] . 

2 The relationship between the scalar integrals in this paper and those in ref. fig ]: h(s) — /2(0) + 2— Wi(s), 
h(s) = -W2(s)/2s and h(s,t, m 2 H ) = -W3(s,u,t,m 2 H ) / st. 
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A. 2 Higher dimension integrals 

The higher dimension integrals are defined in eq. (|A.1| ) with D set to the appropriate value. 
The relationship of these integrals to the usual four-dimensional ones has been extensively 
discussed in ref. [26|| . For n < 6, (6 — 2e)-dimensional integrals can be written in terms of 
(4 — 2e)-dimensional integrals via the integral recursion relation 



rD=6-2e 



(n — 5 + 2e)c 
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D=4-2e 
n 



/ j t ! i n-l 



(i),D=4-2e 



(A.8) 



where q = J2j=i , Co = J27=i c « > an d the matrix 5^ is 



m2 _ > witn P« = , and pij = pji = ki + k i+ i H h kj-i for i < j . (A.9) 



Higher dimension integrals arise naturally when performing the calculations described in 
this paper. In a typical primitive amplitude discussed in the text we obtain integrals of the 
form 



(2tt) 4 (2tt)- 2 ' (p 2 - /i 2 - m 2 ) ■ ■ ■ ((p - YJCi h) 2 - V 2 ~ m 2 ) ' 

(A.10) 

where we have explicitly broken the (4 — 2e)-dimensional momentum into a four-dimensional 
part, p, and a (— 2e)-dimensional part, fi. 

As noted by Mahlon 0, using the definition ( |A.1|) , we can write 



rr>=4-2er..2 



rt 



IS 



-el 



D=6-2e 



(A.11) 



Note that although the loop momentum is shifted to higher dimension, the external momenta 
remain in four dimensions. 



A. 3 Feynman parameter reduction 



As discussed in the text, tensor integrals may be evaluated using Feynman parameters. We 
make use of the following formulae in the text. 

After Feynman parameterization and the normal shift of momentum we make use of the 
following, 

d 4 q d~ 2e fi 



q a i . . . q a 2r+i 







(27r) 4 (27r)-^ ( q 2 - ^2 _ S .. a . a .^n 

where is defined in eq. ( |A.9|) . 

After integrating out the loop momenta we obtain integrals of the form 

/Ofc) 



(A.12) 



I°[f(a k )] = T(n - D/2) C d n a k 5(1 - £ r a r ) 

Jo 



n-D/2 



(A.13) 
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An integral reduction formula that is quite useful is [^6| 



, t U),d 

■ij 1 n-l 



+ ^ 1° 

C » 



(A.14) 



where = S^ 1 — ^ , and q and c are defined in appendix |A.2| . Another useful integral 
reduction formula for I® [a^j] can also be found in 



B Rules for /i's 

To begin, we have a (4 — 2e)-dimensional vector, Q a , which we can express as a sum of a 
four dimensional piece q a and a (— 2e)-dimensional fi a ; Q a = q a + fi a . (We have used the 
convention that upper-case momenta are (4 — 2e)-dimensional and lower-case momenta are 



four-dimensional.) Here we summarize the rules for dealing with /z's given in ref. [10|. 

We abide by the conventional Dirac algebra, {7", 7^} = 2r] al3 and 7^7° = 7°7 Q , 
where a and (3 are (4 — 2e)-dimensional Lorentz indices. The metric is given by rf 1 ^ = 
diag(+, — , .). It follows that Q ■ Q = q 2 — /i 2 (the negative sign preceding the /z 2 is 

from the metric). 

Thus (/l anti-commutes with four-dimensional 7 matrices, = 0. So if we use the 

conventions of 't Hooft and Veltman p8| , and adopt the arbitrary dimension definition, 
75 = i7°7 1 7 2 7 3 , then {^,75} = and ^,75 = 0. That is ^ freely commutes with 75, since 
75 is constructed as a product of four- dimensional Dirac matrices. Hence, u± </[ = <j[ ^=f an d 
w ±j^ = lf LU} ±, where u± = |(1 ± 75) is the helicity projection operator. In calculations it is 
important to first commute the /A's together before evaluating the remaining four-dimensional 
pieces. 

In the calculation presented in section |] the momentum of the sewn fermion lines is 
(4 — 2e)-dimensional. We will borrow the bra and ket symbols to represent the spinors, but 
as helicity is not a good quantum number we shall not label them with a helicity. Since we 
always sum over all states across the cut, there is no need to define a (4 — 2e) helicity notion. 

In sewing the (4 — 2e)-spinors together we implicitly sum over the spin degrees of freedom, 
that is 

\Q)(Q\ =Q + m = $ + ifi + m, and | - Q){-Q\ = -0 + m = -f - + m . (B.l) 
This notation glosses over the distinction of spinors and anti-spinors. However, as discussed 



in ref. [10], this may introduce an overall sign which can be readily fixed by hand. Al- 
ternatively, one may easily use both particle and anti-particle spinors to obtain identical 
results. 
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